The interplay between real and pseudo magnetic field in graphene with strain 
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We investigate electric and magnetic properties of graphene with rotationally symmetric strain. 
The strain generates large pseudo magnetic field with alternating sign in space, which forms strongly 
confined quantum dot connected to six chiral channels. The orbital magnetism, degeneracy, and 
channel opening can be understood from the interplay between the real and pseudo magnetic field 
which have different parities under time reversal and mirror reflection. While the orbital magnetic 
response of the confined state is diamagnetic, it can be paramagnetic if there is an accidental 
degeneracy with opposite mirror reflection parity. 
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PACS numbers: 73.22.Pr, 71.23.-k, 71.55.Jv, 75.30.Cr 

The recent successful preparation of one-atom layer of 
carbons, graphene [H, Q, has provided the opportunity 
of theoretical and experimental research of relativistic 
physics in nanoelectronics. While quantum dots which 
confine quasi-particles in graphene are basic building 
blocks for its nanoelectronic application, the confinement 
turns out to be nontrivial. It is because, in graphene 
where the quasi-particles are described by massless Dirac 
fermions, they can penetrate large and wide electrostatic 
barriers due to the effect of Klein tunneling In prin- 
ciples, graphene dots can be realized by a spatially in- 
homogeneous magnetic field [4-6], but the required mag- 
netic field for the confinement, however, is unreason- 
ably strong compared to usual electronics application. 
Recently, strain engineering of graphene [7J, 18| attracted 
great attention as an alternative tool for graphene elec- 
tronics because the strain induces strong pseudo mag- 
netic field which guides electrons. Thus, for the strained 
graphene to work successfully in combination with exist- 
ing technologies, it is now important to understand the 
physical properties of the pseudo magnetic field. 

In this Letter, we investigate the relative contribution 
of real and pseudo magnetic field to the electric and mag- 
netic properties of the graphene. We show that reason- 
able size of strain can generate strong pseudo magnetic 
field to form a graphene quantum dot with six chiral 
channels. It will be demonstrated that the different sym- 
metry of real and pseudo magnetic fields give rises to rich 
properties of channel opening and orbital magnetism. 
The pseudo magnetic field appears since the variation 
of hopping energies by elastic strains enters the Dirac 



equation [91-4141] . While the strong confinement is due 



to the fact that the pseudo-magnetic field is very strong 
(~ 10T), the six chiral channels are due to the topol- 
ogy of the pseudo magnetic field, where charged particles 
propagate along the zero-field line. As we will show here, 
the real and pseudo magnetic field have different pari- 
ties under the symmetry operation such as time reversal 
and mirror reflection. From the symmetry arguments, we 




FIG. 1: The vertical displacement of our considered systems 
in the text (a) and the pseudo magnetic fields (b). The radius 
is _R=100nm and the height is /io=20nm. 



prove that while the real magnetic field breaks the time 
reversal symmetry in its Hamiltonian, it does not lift the 
valley degeneracy. We will demonstrate our theory by 
showing orbital diamagnetism of the confined state. It 
will be shown that paramagnetic response is also allowed 
when an partially open state with opposite parity be- 
comes degenerate with the confined state. 

Let us consider graphene where mechanical deforma- 
tion is allowed in a restricted disk shape. This can 
be realized by a circular hole made in substrate be- 
low the graphene sheet, and the deformation is induced 
through external force. In experiments, circularly sym- 
metric strain fields can be applied by an AFM tip or by 
a homogeneous gas pressure acting on graphene below 
the substrate [Tij. When strain is induced by homoge- 
neous load, the optimized vertical displacement h(r)is 
given by [la] 
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where /o is the force per unit area acting on surface, D 
is the bending rigidity, ho = /o-R 4 /(4 3 -D) is the vertical 
displacement at center, and R is the radius of the region 
where the deformation is allowed. The in-plane relax- 
ation of the carbon atoms u x , u y can be calculated by 
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pseudo gauge field A ps (r) is written as 
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FIG. 2: The eigenenergies as a function of the system size L. 
We set i?=100nm and fto=20nm. Note that at certain energy 
(marked by an arrow), the eigenenergy becomes insensitive to 
the system size L . 



minimizing the elastic free energy for the given vertical 
displacement h(r), 



F = / dxdy 



A 



(2) 



where k is the bending rigidity, A and 3A) are Lame 
coefficients, and My is the strain tensor. Here, the strain 
tensor u i3 -(r) is related to the displacement fields via 
u xx = du x /dx + ^(dh/dx) 2 , u yy = du y /dy + \{dh/dy) 2 , 
and u xy = \(du x / dy + du y / dx) + h(dh/dx)(dh/dy). 

We consider spinless fermions in a graphenc lat- 
tice. The spinless quasi-particle in the graphene can 
be described by 4-component wavefunction \& T = 

(^k>^-k) = {i>A,K,^B,K^A,-K,i>B,-K)- These are 
the electron wavefunctions near two inequivalent points 
±K in hexagonal Brillouin zone in the two crystalline 
sublattices A and B. We ignore the valley mixing by the 
strain based on the results of a tight-binding calculation 
showing that the two inequivalent valleys are not coupled 
under uniaxial deformations up to 20% (l7j . The main 
effect of the strain field on the electrons is to modify the 
energy for electron hopping between the nearest neighbor 
atoms. The modified quasi-particle energy by the strain 
is well described by introducing the pseudo gauge field 
Ap S (r) in the massless Dirac equation, 
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(3) 



where — e is the electric charge, vp ps 9 x 10 5 m/s is the 
Fermi velocity, p = —iK(d x ,d y ), and a = (a x ,a y ) are 
Pauli matrices acting in the sublattice space. Here, we 
take two K points in opposite direction along x-axis. The 
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where t w 2.8eV is the electron hopping energy between 
the nearest n orbitals, and /3 ~ 2 — 3 is the dimension- 
less coupling parameter for the lattice deformation B 
Hit ]. In FigJTJ we plot the pseudo magnetic field 
V x A ps . This is the pseudo magnetic 
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B ps k = 



field experienced by the particle in K valley, and for the 
particle in —if valley, the sign of the magnetic field is 
opposite. 

To investigate the confinement of the quasi-particles 
in graphene, we compute the eigenenergies of the sys- 
tem as the total graphene size used in the calculation 
L increases and observe whether the energy is insensi- 
tive to the system size. The graphene in the calcula- 
tion is a disc defined as r < L and the graphene is de- 
formed only in the region of r < R. We impose the 
boundary condition ^^(r = L) = ^s(r = L) = 0. 
While the wavefunctions which extend over the total sys- 
tem are sensitive to the boundary conditions, localized 
states in the deformed are not affected by the bound- 
ary condition. We obtain the eigenvalues of Hamilto- 
nian using basis functions ipA,B 



Y, c tf^riu(r), where 



0n,;( r ) — Ji(^jr-r)e lW . Here, a n i is the n-th zero of the 
Bessel function J; (a;) of order /, and we use the indices 
/ = 0,±1,±2,---, n = l,2,3,---. 

As shown in Fig. 2, we find that at certain energies, 
there exist eigenstates of which eigenenergies are insen- 
sitive to the system size L. These are the localized 
states induced by the deformation of the graphene. Com- 
pared to the mid-gap states in a ripple array studied in 
Ref. Q which shows weak size dependence in logarithmic 
scale, the eigenenergies of the localized states show al- 
most no size-dependence. From an analytic calculation, 
the length and energy scales for the localized state can be 
obtained by bringing the asymptotics of differential equa- 
tion for the wavefunction to dimensionless form. This is 
done by rescaling r — > rpf, E —¥ EqE, with 
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E Q = hv F /r 



(5) 



The scale ro thus plays a role of the localization length 
and can be estimated as (R 4 a/h 2 ) 1 ^ 3 , where a is the 
lattice constant. Since the localization length must be 
shorter than the hole radius, otherwise the Dirac equa- 
tion with effective magnetic field can not apply, the lo- 
calization length expressed in the above equation is valid 
only for strong enough load, ho 3> %/ Ra. The scale Eq is 
associated with the depth of the potential well and the 
energy of the localized energy levels. 

It proves useful to consider the symmetries to under- 
stand the energy spectra of quantum systems. The time 
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reversal symmetry is not broken by the strain if we con- 
sider the problem with both of the valleys. A time rever- 
sal operation defined by 



T 



(6) 



satisfies THT^ 1 = H, TiT' 1 = -iyM. 

Here K, is the complex conjugate operator and 
(t x , T y ,T z ) are the Pauli matrices acting in the the valley 
space. Note that, the Kramers degeneracy is not relevant 
here since the original system has orthogonal symmetry 

r 2 = i. 

The quasi-particle can stay in a given valley provided 
there is no short range scattering ( e.g. lattice defects). 
In this restricted case, the single valley time reversal sym- 
metry is broken by the pseudo magnetic field. The single 
valley time reversal operator for 1/2 spin S — —ia y JC 
does not commute with the Hamiltonian in Eq.Q. The 
Hamiltonian is symmetric under the symplectic time re- 
versal transformation only in the absence of the strain 



(S 2 



-1); SH{A ps = 0)S- 1 = -H(A ps = 0). The 



relevant Kramers degeneracy here is lifted by the pseudo 
magnetic field A ps ^ 0. 

The Hamiltonian in Eq.@ is also symmetric under a 
mirror reflection 



M = 
MHM~ 



H, 



(7) 
(8) 



where k x acts as i - > x, y — > —y. The valley in- 
dex remains same under the mirror reflection but in- 
evitably changes lattice index. The spatial symmetry 
(\iPa,k{x, y)\ 2 = \ip B ,K(x,-y)\ 2 ) of the probability den- 
sity of the localized state shown in Fig. 3 (a) and (b) 
reflects the AJ-symmetry in the Hamiltonian. 

To understand the electronic structure of the strained 
graphene, we rewrite the Hamiltonian in Eq. (|3j) for a 
given valley K and eigenenegy E, 



v 2 F (p + eA ps ) 2 + hev 2 F (W x A ps ) z cr 2 ^ K = E 2 V K .(9) 



The first term in the left side of Eq.® comes from the 
kinetic energy and the second term is due to the pseudo 
Zeeman coupling. The pseudo magnetic field is strongest 
at six points forming a hexagon (Fig. 1) where local Lan- 
dau levels might be formed. For a given pseudo spin and 
valley, one can see maximum probability density around 
only three points (Fig. 3(a) and (b)). This is because of 
a necessary condition for the stable confinement on the 
pseudo Zeeman coupling; 



(V X A ps ) 2 cr z < 0. 



(10) 



In this case, it costs higher energy when a quasi particle 
to go out to weaker field region. The triangular (instead 
of hexagonal) shape of the wavefunction for a given lat- 
tice and valley is due to the selective stabilization by the 
pseudo Zeeman coupling. 




FIG. 3: (a) Probability densities of the localized wave- 
functions \tp A ,K(r)\ 2 and (b)|i/> s , K (r)| 2 . We set i?=100nm, 
fto=20nm. (c) Classical closed orbits of a charged particle 
in the inhomogeneous magnetic field given by B pa which re- 
semble the probability density of the localized states shown 
in (a) and (b). (d) Classical paths describing three outgoing 
channels. 



The confinement of channeling in the strained 
graphene can be visualized by investigating classical tra- 
jectories of the charged particles in the pseudo magnetic 
field (Fig. 3 (c) and (d)). Among the periodic orbits 
around the pseudo magnetic field maxima, we find clover- 
shaped orbits which resemble the localized wavefunctions 
in (Fig. 3 (c)). These closed orbits are very unstable 
against small perturbation. In quantum mechanics, the 
clover-shape motion might be responsible for quantum 
transition between the sites of the local density maxima. 
We also find out-going trajectories (Fig. 3 (d)) for differ- 
ent initial velocities. A charged particle can propagate 
along the line where magnetic field changes sign, which 
is, so called, snake orbit\2ty. Due to the symmetry of the 
pseudo magnetic field, there are incoming trajectories in 
60 degree rotated angle from those of outgoing trajec- 
tories. Quantum mechanically, for given components of 
lattice and valley, the graphene quantum dot is connected 
to three incoming and three outgoing chiral channels. 

The opening of the channels manifest in the energy 
spectra. The localized energy level undergoes crossing 
and avoided crossing as the graphene size L changes. 
( see Fig. 2 and Fig. 4 (a) for mode detail.) When 
the channels open, the eigenenergy of the confined state 
is affected by the graphene size, so has avoided cross- 
ing. Meanwhile when the two level cross each other, the 
channels remain closed and the eigenenergy of the local- 
ized state is insensitive graphene size. In the process of 
avoided crossing, the confined state undergoes transition 
to an outer state and a new outer state becomes local- 
ized. Since any parity does not change in this continuous 
process, the confined state transit only to the state with 
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same parity. 

Let us consider the response of the strained graphene 
to real magnetic field. When real magnetic field, B rc = 
V x A ro , is applied, the minimal coupling of the elec- 
tromagnetic gauge field A re is done by replacing p with 
p + eA r eai in Eq. (|3]) . The Hamiltonian becomes H + H' 
where 

H' = v F e (jj _°^) ■ A ro , (11) 

and A ro = ^(—y,x)B is the gauge field for the real mag- 
netic field in z-direction. We want to address here that 
the application of real magnetic field breaks the time re- 
versal symmetry of the Hamiltonian, but it does not lift 
the valley degeneracy. This can be proved by showing 
< ^±k\^-\^±k >= 0. The proof comes from the fact 
the eigenstates have either even or odd parity of the mir- 
ror reflection symmetry in Eq.©, Ai 2 = 1: 

dH ' 

< *K|-^g|*if >=< ^K\(~ycr x + xa y )\y K > 

= < y K \M(+ya x - xa y )M\^ K > 

= - <y K \(-ya x + xa y )\y K > . (12) 

Since < ^k\^-\^k > is equal to its own minus value, it 
must be zero. The leading magnetic field dependence of 
the eigenenergy in the presence of magnetic field is not 
linear but quadratic, oc B 2 . It comes from the kinetic 
energy and its sign is positive. Therefore, the orbital 
magnetization of the strain-induced quantum dot at zero 
temperature, is diamagnetic ( —dE/dB < ) and pro- 
portional to the applied magnetic field strength (See Fig. 
4 (b)). 

In contrast to the diamagnetic response of the confined 
state, the orbital magnetic response can be paramagnetic 
( —dE/dB > ) when there are level-crossings. Near 
the region of the level crossings, there are two energy 
levels with opposite parity of Ai. One of the states is a 
localized state and the other is a partially opened state 
(not shown). The accidental degeneracy which happened 
here can be lifted by applying real magnetic field because 
H 1 is odd under the mirror reflection, M.H'AA = —H'. 
Then energy splitting arises, proportional to the real 
magnetic field strength, which contributes to paramag- 
netic response. 

In conclusion, we have shown that the rotationally 
symmetric strain in graphene can be considered a quan- 
tum dot with spatially separated six chiral channels. The 
chiral channels exist along the line where the pseudo mag- 
netic field changes sign. The real and pseudo magnetic 
field has different symmetry under a mirror reflection, 
which makes the orbital magnetism be diamagnetic or 
paramagnetic depending on the degeneracy. The orbital 
magnetic response of the confined state is diamagnetic 
due to its kinetic energy. When there is an degeneracy 
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FIG. 4: (a) A detailed energy spectra as a function of the 
system size L which shows crossing and avoided-crossing. (b) 
The orbital magnetic response, — §§ , of localized states. 
Mostly the magnetic response is diamagnetic — |j| < but 
the response for a level crossing point (indicated by dotted 
line) shows paramagnetic response (positive value). We set 
i?=100nm and /io=20nm. 

with opposite mirror reflection parity, the orbital mag- 
netism can be paramagnetic. 

Quite recently, we became aware of a work on dynamics 
of electrons in strain- induced pseudo magnetic fields plj]. 
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